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Abstract 

A system consisting of a chaotic (billiard-like) oscillator coupled to a linear wave equation in the three- 
dimensional space is considered. It is shown that the chaotic behaviour of the oscillator can cause the 
transfer of energy from a monochromatic wave to the oscillator, whose energy can grow without bound. 



1 Setting the problem and result 



The system we consider (a linear wave coupled to an oscillator) is formally defined by the Hamil- 
tonian 

h = \(p 2 v+ pi) + v & z ) + ( Mv, z) f u(x, t) d 3 x + \ f (u 2 + (V x u) 2 ) d 3 x, (1) 

2 J\\*\\<1 2 J 

where u(x, t), the massless Klein-Gordon field, is a scalar function on R 3 xt 1 , and (y, z) £ M 2 are 
coordinates in the configuration space of the oscillator ((p y = y,p z = z) £ R 2 are the corresponding 
momenta). The smooth potential V(y, z) is bounded from below and tends to infinity as ||y, z\\ —¥ 
oo. To avoid technicalities, we assume that V equals to infinity outside a bounded domain in 
the (y, z) plane. The interaction coefficient k(y, z) is smooth and bounded along with the first 
derivatives, and e is small. 

The corresponding equations of motion are 

Utt-Au = -ek(y(t) ) z(t))Z(x.) ) (2) 
where £(x) = | J JJ X JJ ~ j (the characteristic function of the unit ball in the x-space), and 

V + i~ (V(y, z)) = -ek' y [ u(x, t) d 3 x, 

9V 7||x|1 ^ (3) 



z+^{V{y,z)) = -ek' z [ «(x,()d 3 x. 
oz J\\x\\<l 



We may think of the oscillator as being located in the unit ball in the x-space and emitting/receiving 
the wave w(x, t). It can be shown (similar to [TJ HJ |3]) that if the initial energy of the wave is finite 
(i.e. J [u 2 + (V x u) 2 ) d 3 x < oo at t = 0), then the field u tends to a constant and, for a typical 
choice of the interaction coefficient k(y, z)0 the energy of the oscillator decreases, and it comes to 



1 For a general choice of k(y, z) the state of oscillator may approach an invariant set on which the function k(y, z) 
stays constant, see more in [HE] 
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a rest at some stationary point of the potential function as t — >• +00, i.e. the energy flows from 
the oscillator to the field and is carried away to infinity. In this paper we show that if the wave 
has infinite energy, then an opposite process may take place. Namely, when the oscillator operates 
in a chaotic regime, the energy can be pumped from the field to the oscillator, and the oscillator's 
kinetic energy can increase up to any given value. 
The field equation ([2]) can be explicitly resolved: 

«(x, t) = uo(x,t) + ±J mt - "y - ! ' Sil)) e(x - s) d 3 s, (4) 

where u is the solution to the homogeneous wave equation. We take as u a a monochromatic 
standing wave of infinite energy: 



Uq(x, t) = Ksm(u)t) J a k cos(k • x) S(ui — ||k||)e? 3 k 



(5) 



(the fact that it is a standing wave is not very important; the fact that the frequency spectrum 
is discrete and finite is used in an essential way). By placing (j4]) into ([3]) we obtain a system of 
delayed differential equations 
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y+^(V(y,z)) = - A^k' y {y,z)sm{ut)-e 2 k' y J k(z(t - s),y(t - s))P(s) ds 



2 



{V(y, z)) = - A u ek' z (y, z) sm(ut) - e 2 k' z / 



(6) 



where 



A u =kJ a k £(x) cos(k ■ x ) 5(w - j|kj| )d 3 kd 3 x, (7) 

and V(s) = — (16 — 12s + s 3 )s (we denote here s = ||s|| ). This is a Hamiltonian system, subject to 

a conservative periodic perturbation of order e and a dissipativ^ correction of order e 2 . Namely, if 
we drop the 0(e 2 )-terms (those including the delay), system (|6|) will be defined by a time-dependent 
Hamiltonian function 

H = l - (p 2 y + pi) + V(y, z) + eA^k{ yi z) sin(^). (8) 

In general, the systems with a time-dependent Hamiltonian do not preserve energy, and the periodic 
forcing may, in fact, lead to an unbounded growth of energy. Importantly, as the kinetic energy 
grows, the forcing becomes effectively slow. It is known [H [5J [6j [7] that if a Hamiltonian system 
behaves chaotically at all sufficiently large energies, then adding a slow periodic forcing creates 
orbits of unbounded energy growth. Similar to [7] , we will show that in the system of type ([8]) this 
growth is linear in time; namely, one may construct orbits for which the energy gain per period of 
the force is bounded from below by a non-zero constant of order e. Using the fact that the delay 
terms which can lead to a dissipation of energy are of order e 2 , i.e. they are much smaller than 
the energy gain, we show below that the full system ((6]) has, for sufficiently small e, solutions for 
which the energy grows up to any given value, linearly in time. 

Our method does not allow us to show the existence of orbits for which the energy tends 
to infinity (the larger the energy value we want to achieve, the smaller value of e we have to 



2 We are not showing here that the delay term always leads to a dissipation of energy. In order to get some 
insight, one can check that in the linear case (i.e. quadratic potential V and linear coefficient of interaction k) the 
addition of this term shifts the spectrum to the left of the imaginary axis, i.e. creates dissipation (however, as we 
consider bounded functions k and the potential V that is infinite outside a bounded domain, the linear case is not 
the subject of this paper). 
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take). Similarly, our approach is not applicable to the case of non-zero mass in the wave equation. 
However, we conjecture that such orbits do exist in system like they do in the time-dependent 
Hamiltonian system defined by (JSJ (see [7]). We also think that our approach (based on a reduction 
to an invariant manifold) is interesting in its own right. We stress that the mechanism of the 
sufficiently fast (linear) energy growth we obtain here is based on the chaotic behaviour of the 
oscillator. Would the oscillator be integrable, the existence of adiabatic invariants |9] would 
impede a possible energy growth in the shortened system (jHJl, and it seems plausible that the 
0(e 2 )-dissipation due to delayed terms should arrest the energy growth completely in this case. 
So, the picture we propose is the following: while an integrable oscillator that interacts with a 
monochromatic wave (in R 3 ) should lose energy, a chaotic oscillator may take the energy from the 
wave. Note that the reported effect cannot be immediately interpreted as a resonant phenomenon: 
since the potential is infinite outside a bounded domain, the region of allowed motions (the Hill's 
region) in the configuration plane (y, z) is always bounded, and since the velocities p y and p z tend 
to infinity as the energy grows, it follows that the characteristic return times tend to zero and 
become much smaller than the period of external force. 

In order to formulate the result precisely, we need to define what do we mean by the "chaotic 
oscillator" . Consider a system of ODE's defined by the Hamiltonian 

H=±(p 2 v +p 2 z )+V(y,z) (9) 

(the Hamitlonian (j8]) at e = 0). As we mentioned, we assume that V is infinite outside a certain 
bounded region D. Following [7], assume that for each sufficiently large h, the system defined 
by (0) has, in the energy level H = h, a pair of hyperbolic periodic orbits L a and Lf, such that 
the unstable manifold of L a has an orbit r a f, of transverse intersection with the stable manifold 
of Lb, and the unstable manifold of Lb has an orbit Tb a of transverse intersection with the stable 
manifold of L a . Moreover, we assume that the orbits L a , Lb, ^ab and Tb a depend continuously 
on h. When these conditions are fulfilled we call system ([9]) a chaotic oscillator. It is well known 
[TU1 ITT) that the existence of the transverse heteroclinic cycle implies a chaotic behaviour indeed. 
Namely, fix a value of h, and take a sufficiently small neighbourhood U of the heteroclinic cycle 
L a U Lb U r a 6 U Tba in the energy level H = h. Then the set of all orbits that stay in U (~l {H = h} 
is in one-to-one correspondence with arbitrary sequences of the symbols a and b (one round made 
by the orbit near L a is coded by a, and a round near Lb is coded by b). 

The basic example of chaotic oscillators in our setting is given by systems with billiard-like 
potentials. Namely, let D C IR 2 be the bounded domain such that V(y, z) is finite inside D and 
infinite outside of D. Let the boundary of D consist of a finite number of smooth arcs, Si, . . . , S n , 
joined at corner points. Let C r -smooth functions Qi(x,y), . . . ,Q n (x,y) be such that, for each 
j = 1, ... ,n, the function Qj is defined in a neighbourhood of the boundary arc Sj, the arc Sj is 
a level line of Qj (i.e. Qj{x,y)\( x ,y)£Sj = const), and VQj ^ in the neighborhood of Sj. We 
will call the potential V that equals to infinity outside of D a billiard-like potential if there exist 
C r -smooth, strictly monotonic functions Wi(Q), . . . , W n (Q) such that for each j = 1, . . . , n the 
potential V in a small neighbourhood of the arc Sj is given by 

V(x > y) = W j (Qj(x,y)) (10) 

(we do not include the corner points into the arcs Sj, i.e. they are open intervals, so their small 
open neighbourhoods do not need to contain the corner points, hence even if two arcs join at a 
corner point, their small neighbourhoods where ((TO)) holds do not need to intersect, i.e. no relation 
between the corresponding functions Wj arises). By scaling the momenta p VtZ to \fh, system (j9)) 
on the energy level H — h transforms into the system 

H=\{pI+P 2 z ) + \v{ V ,z) (11) 



3 



on the energy level H — 1. One can easily check that if V is a billiard-like potential (i.e. it is defined 
by (JTUJ) near the boundary arcs), then the family of Hamiltonians (fTTj) (with hr 1 being a small 
parameter) satisfies conditions of |14j which guarantee that the flow defined by such Hamiltonian 
is an approximation, at h large enough, to the billiard flow in D. Namely, billiard in D is the 
mechanical dynamical system which describes the following motion of a point mass |18j : the 
particle moves inertially (with a constant velocity) inside the region D in a plane until it hits the 
boundary of D, then the particle is reflected according to the elastic reflection law, the angle of 
reflection is the angle of incidence, and so on. As it follows from [14) . the time-shift maps by the 
billiard flow in D and by the Hamiltonian flow defined by (fTTj) with a billiard-like potential V and 
h large enough are C r -close outside the set of singular billiard orbits (i.e. such billiards orbits 
which enter a corner point or are tangent to the boundary of the billiard domain D). Therefore, 
if the billiard in D has a hcteroclinic cycle of two hyperbolic periodic orbits L a and Lj, and two 
transverse heteroclinic orbits T a i, and Tt, a and neither of these orbits is singular, then system 
in the energy level H = 1 (so system (0) in the energy level H = h) also has such heteroclinic 
cycle for all h large enough, i.e. system ([9]) with a billiard-like potential V is a chaotic oscillator 
according to our definition. In particular, the oscillator defined by a billiard-like potential is chaotic 
when the underlying billiard is dispersive: all the boundary arcs are concave and meet each other 
at non-zero angles at the corner points [I5j [HI H3 HB1 E] ■ 

Let L a and Lf, be the two non-singular hyperbolic periodic orbits of the billiard in D that are 
connected by the transverse non-singular heteroclinics r a f, and T{, a - Fix the speed of the particle in 
the billiard to be equal to 1, and let T c (where c = a,b) be the period of L c , and (y c (t), z c(t))\tG[o.T c ] 
be the equation of the orbit L c . Denote 

f f T ' 

Vc = 7 F / k{y c {t),z c {t))dt. (12) 

J c Jo 

Assume 

v a v b . (13) 

Theorem 1.1. Consider a system (0), $3$) that describes a linear massless scalar field interacting 
with a chaotic oscillator, with a billiard-like potential V and the interaction coefficient fc(y, z) 
bounded with first derivatives. Let condition (T^) hold. Then there exists ho such that for each ui 
for which A M ^ (see ^) and for any hi > ho there exists eo(h\) > such that for all e 6 (0, eo) 
the system has a solution with the wave component u given by |7p,([I|j and the oscillator component 
{y{t), z(t),p y (t),p z (t)) reaching from the energy ho at t = to the energy h± at some finite t (the 
energy of the oscillator is given by 

The proof of the theorem occupies the rest of the paper. In Section [5] we prove an invariant 
manifold theorem which allows us to reduce the system of delayed differential equations © to 
a four-dimensional non-autonomous system of ODE's, which is 0(e 2 )-close to the Hamiltonian 
system given by ©. This reduction is possible on any bounded set of values of (y, z,p y ,p z ) for 
sufficiently small e. Note that the maximal value of e for which we can guarantee the reduction 
increases as p y and p z grow, which is one of the reasons why we show only bounded energy growth 
in Theorem 11.11 In Section [3] we apply the construction of [7] (modified for the non-Hamiltonian 
case) to the reduced system, and finish the proof. 

2 Invariant Manifold Theorem 

In this section we reduce the infinite-dimensional system of differential equations with delay to a 
finite-dimensional system of ordinary differential equations. Consider the system 

X(t) = F(X(t))-S [ G(X(t - s), s) As, (14) 
Jo 
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where X belongs to an n-dimcnsional smooth manifold M, and the functions F and G arc C r - 
smooth. Let F be such that the differential equation 



X = F(X) (15) 

has, for any initial condition, a solution defined for all t € (—00, +00) (for example, this differential 
equation is Hamiltonian with compact energy levels). System ([6]) can be represented in this form 
with S = e 2 . Namely, we introduce a variable 6 £ S 1 and rewrite ([6]) as 



V=Py, z=Pz, v = u 
d 



Py = ~jr (V(y, z)) - A^ek'yiy, z) sin(0) - e 2 k' y I k(z(t - s),y(t - a))P(a) ds, 

dy v Jo (i 6 ) 

8 



2 



(i 



2 

Pz = -^-{V{y,z))-A u ek' z (y,z)sm(6)-e 2 k z \ k{z(t - s),y{t - s))V{s) ds, 







so X = (y, z,p y ,p z ,8 = tot) 6 I 4 x S 1 here. 

It is well known that given any continuous function X : [— r, 0] — > R n there exists a unique 
solution X(t) of equation p^|) such that X(t) = X(t) at t £ [— r, 0]. One can therefore view the 
evolution defined by equation (| 14[) as a semiflow in the space C of continuous functions that act 
from [— r, 0] to R n : the time-s map of the semiflow takes the initial condition X(s)|[_ t ] to the 
segment X(s)|[ t _ T) t] of the corresponding solution. A smooth function \i : R n x [0,r] — > R n defines 
a map R n — > C by the rule X(s) = n(x, —s); the graph of such map is an invariant manifold for 
the semiflow defined by (|14l) if the solution with the initial condition fj,(X(0), — s)| se r_ Ti0 ] satisfies 

f x(X(t),s)=X(t-s) (17) 

for all t > and s e [0, r]. When such manifold exists, the restriction of system (|T4|) onto it is a 
system of ordinary differential equations 

X = F{X) - S [ G(n(X, s), s) ds. (18) 
Jo 

In other words, the existence of a smooth function fi which satisfies the invariance condition (fTT)) 
implies the existence of an n-parameter family of solutions (parameterised by -X"(0)) to the delayed 
differential equation (TH)) which also solve the ordinary differential equation pi 



Lemma 2.1. Given any compact subset K 0/ the X-space M, for all sufficiently small 5 there 
exists a C r ~ x -smooth function [i : K x [0, r] — > R n such that for any X(0) £ K the solution X(t) 
of Jl^[ ) which starts with the initial condition n(X(0), — s)| se r_ r , ] satisfies the invariance condition 
Ji7| ) for the interval oft values for which the solution stays in K. 

Proof. By the Cauchy-Picard-Lindclof theorem }19) . given any smooth functions fi, F and G, the 
ordinary differential equation (|18|) generates a uniquely defined (on any given finite time interval) 
solution for any initial condition from K if 6 is small enough. This solution depends smoothly 
on any parameter on which the system depends smoothly. In particular, the solution depends 
smoothly on the function /i. Thus, we will show below that the solution is a C^-function of [i 
which is considered as an element of the space of C"" -1 -smooth functions. 

Given /1 and S, take any X(t) £ K, consider its backward orbit by equation (fTSJ), and let 
X(t — s) be the point on this orbit which corresponds to the (backward) shift to time s. Denote 
as 4>(/j,,d) the map K x [0,r] — > R n which (for given [i and 6) sends X(t) £ K and s £ [0,r] to 
X(t — s), i.e. 4>{pL, S) is the backward flow of the ordinary differential equation (fT51) . The flow has 
the same smoothness as the equation, so if fi £ C r ~ 1 , then 4>(/j,, 6) is C r ~ 1 with respect to X and 
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s. At 6 = the flow <j>(n, 6) is generated by equation (|T5|) . so it is independent of //. We will show 
in a moment that <p depends C^-smoothly on /i and S. As 4>(n, 0) does not depend on /i, it follows 

that the Frcchct derivative — vanishes at 5 = 0. Hence, by the Implicit Function Theorem, the 

a/j, 

equation 

H = 4>(ji,S), (19) 

has, for every small S, a unique solution /i € G r_1 . If we plug this particular yu into the right-hand 
side of (fT8| . then equation (fT9| will exactly mean that condition (fT7|) is satisfied by the solutions 
of (fl8|) . Thus, each of these solutions will also solve the original delayed equation (fT4| . which is 
the statement of the lemma. 

Thus, to finish the proof, it remains to show the smooth dependence of the flow of (fT5)) on fi. It 
is well known that the solutions of ordinary differential equations depend smoothly on the function 
in the right-hand side of the equation. So we are left to show that the so-called Nemytsky operator 
(or substitution operator) Af which takes the function /j,(X, s) to the function G(/i(A, s), s) is of 
class G 1 on the space of G r_1 -functions /x, provided G is C r as a function of (X, s). In order to do 
this, it is enough to check that the Frechet derivative of N at a given function \x is the operator of 

QQ 

multiplication to -x— (/i(A, s), s). For that, one needs to check that 
oX 

+ A/j) - G(fx) - G'(fx) ■ A/iH^ = oGIA/iU^) 

(we suppress, notationally, the dependence of s, so G' denotes here the derivative of G with respect 
to its first argument). This relation is rewritten as 



[ (CfiH + tAn) -G'in)) dC-A^ 
Jo 



= (||A M || cT ._ 1 ) 

C r-1 



which reduces to the obvious (since G is C r ) claim that G'(/j + £A/x) — ^G' (/i) as Aix^-^^, 
uniformly for all £ G [0, 1]. □ 



Applying this lemma to system (|16j). we find that for all sufficiently small e it has a family of 
solutions which satisfy the system of ordinary differential equations 

d 

y+fy z )) = - A * e K(y> z ) sin ( 61 ) + e2j: ^y^ z ^ y> z > e > e )> 

z + L (V(y, z)) = - A^k^y, z) sm(6) + e 2 F 2 (y, z, y, z, 9, e), (20) 
6 = u), 

where the smooth functions incorporate the delay terms. Moreover, for all sufficiently small 
e, every solution of this system satisfies the original delayed equations (O on the time interval the 
solution stays in a bounded ball in the phase space (we may take this ball as large as we want; 
however, to increase the radius of the ball, we might need to take e smaller). 



3 A partially-hyperbolic set and the energy drift 

As we just have shown, it is enough to establish the existence of the solutions of growing energy in 
system ([20)) . By taking a sufficiently large ho and scaling time to \fh~o~ : we rewrite this system as 



G 



19 e e 2 — 

V + T-n- (V(y,z)) = - A w — fe (y,z)sin(0) + —Ti{y,z, y/h Q y, y/h z, 9, e), 

1 <9 e e 2 

I + t "o (V(y,z)) = - A ul —k' z (y,z)sm(0) + —F 2 (y,z, y/hoy, \/ho~z,6,e), ( 21 ) 
tio oz Hq iiq 



Vho 

At e = the first two equations are independent on 9: 

Yl (22) 

As we mentioned, the assumption that our billiard-like oscillator is chaotic means that this system 
has, at sufficiently large ho, a uniformly- hyperbolic set A in every energy level H > 1 (where the 
rescaled energy H is given by with h = ho). Namely, we take the two non-singular hyperbolic 
periodic orbits L a and Lb of the billiard in D and the two transverse non-singular heteroclinics 
T a b and Tb a that connect them. As we mentioned, by virtue of |14j . the hyperbolic hcteroclinic 
cycle persists in the smooth Hamiltonian approximation (|22[) of the billiard, provided ho is large 
enough. This means that system (|22[) has in every energy level H = h > 1 a pair of hyperbolic 
periodic orbits L a (h) and Lb(h) and the heteroclinic orbits, T a b(h) and Tb a (h) where r a f, is the 
transverse intesection of the unstable manifold of L a with the stable manifold of Lb, and similarly 
for T ba - 



r ab c W u (L a ) n w s {L b ) 7 r ba c W u (Lb) n (L a ), 

and these four orbits are close (at ho large enough) to the corresponding orbits of the billiard in 
the same energy level. 

By [10], the set A(h) of all orbits which stay in a small neighbourhood (in the level set H = h) of 
L a (h) U Lb(h) U T a b(h) U rb a (/i) is a uniformly- hyperbolic set which is in one-to-one correspondence 
with the set of all sequences of a's and 6's. More precisely, in each energy level we take two small 
smooth cross-sections, E a and £&, to L a and Lb respectively ( we assume the cross-sections depend 
on h smoothly). Every orbit from A must intersect E a U £& infinitely many times. The sequence of 
a's and 6's denoting the corresponding cross-sections that the trajectory passes through gives the 
code, {£i}it°^oo (£i € { a ifr})j °f the trajectory. Namely, if Mj is the sequence of points at which 
an orbit from A intersects E a U £&, then M, £ S^. 

The flow induced by Q22[) on any given energy level H = h defines Poicare maps on these 
cross-sections, denoted by 

II CC , : E c /, c€{a,6}. (23) 

These maps are such that for a trajectory of ([22]) that intersects E a U at a point Mj S £ c and 
then at a point Mi+i £ £ c / we have 

M i+1 =n cc /Mi. (24) 

In other words, II aa and Hbb are the Poincare maps near the periodic orbits L a and, respectively, 
Lb, while n a b and IIb a correspond to a passage near the heteroclinics r a b and r& a . The theory 
of billiard-like potentials built in [14] implies that since the billiard orbits L a , Lb, T a b, Tb a are 
non-simgular, the Poincare maps II CC ' for system ([22]) are C-close, at ho large enough, to the 
corresponding Poincare maps defined by the billiard flow. In particular, the hyperbolicity of these 
maps for the billiard flow (which follows from the fact that the periodic orbits L a and Lb are 



7 



hyperbolic and the heteroclinic orbits r a £, and Tb a are transverse) is inherited in system (|22[) for 
all ho large enough. 

In fact, since each periodic orbit is saddle, one can represent the two-dimensional cross-sections 
£ aj 6 as the cross product of two ceratin small intervals 

E = U a x W a , Z b = U b x W b , 

where U ai b £ K correspond to contracting directions and W a ,b S K correspond to expanding 
directions. By virtue of since the heteroclinic orbits r a f, and Tb a are transverse, one can 
write the Poincare maps (|23l) in the so-called cross form (see [12M13] ). Namely, there exist smooth 
functions 

fee : U c x W d -)■ U' c , g cd : U c x W d ->■ VF C , (25) 

(where c, c' G {a, &}) such that a point Mj = (itj, Wi) € E c is mapped to Mi+i = (wj+i, u^i+i) £ S c / 
by the map n cc < , if and only if 

f CC '{ui, w i+ i) = Ui+i, g C c'(ui,Wi + i) = Wi. (26) 

Moreover there exists A > such that 

d(fcc',9a 



d(u, w) 



< A < 1. (27) 



The latter inequality means, essentially, that the Poincare maps II CC < are contracting in the u- 
coordinatc and expanding in the u-coordinate. 
From [7] it follows that the operator 

{K.Wi)}£^oo -> {(ki-iii{ u i-U w i)>9tei+i( u i> w i+l))}i^oo> ( 28 ) 

is a contraction mapping for any code sequence £ = {^lit"^^, which gives us the existence and 
uniqueness of the orbit L^(h) e A(h) with the given code £ (see [TUl HI])- By fixing a code £, the 
family of the orbits with this code parameterised by h, i.e. the manifold 

l<h<hi/h 

for any given /ii > hi, is a normally- hyperbolic invariant manifold of system (|22p . System (|21l) at 
e = is obtained from ([2^)1 simply by adding an equation = ui/ \fh~o~ for the phase 9 £ S 1 . As this 
equation is decoupled from the first two, and the evolution of 6 is non-hyperbolic (it is just a linear 
rotation), the manifold £° x S 1 is the invariant normally-hyperbolic manifold for system (|21j) at 
e = 0. The normally-hyperbolic invariant manifolds are known to persist at small perturbations 
[20] , which implies that (|2Tj) has an invariant manifold (e) close to x S 1 for all small e (one 
might need to take smaller e to make hi larger). Because the normal hypcrbolicity is uniform, these 
manifolds exist for all codes £ for the same range of e values. A formal proof of the existence of 
the invariant manifolds £j(e) is achieved as follows. Note that the Poincare maps H cd arc defined 
for all e sufficiently small and depend smoothly on e, so they can be written in the form 

Ui+i = fcc'(v*i,w i+ i) + 0(e), Wi = g CC '(ui,w l+ i) + 0(e), 

1 (29) 



h i+ i = hi + 0(e), e i+1 =6 l + 0(e + 



/ho' 



where the energy h is defined by ((SJ; the 0(-)-terms are functions of Ui,Wi+i,hi,9i. Now, lemma 
1 of [7] is applied to these maps, which immediately gives the existence of the invariant manifolds 
£{(e) close to x S 1 provided e is small and ho is large enough. 
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By differentiating (|5J), we find that the rate of change of h = H(y, z,p y ,p z , 9) along an orbit of 
(j2"Tj) is given by 

^ =e(LjA w k(y,z)cas0 + O(e)). (30) 

Since the Poincare map (|29|) at c = d corresponds to one round near the billiard periodic orbit L c , 
it follows from (|30[) that given an arbitrarily small i5 one can choose e and h^ 1 small enough such 
that 

h i+1 ~ h Q l >e(A w v c co80-S) (31) 

in the map IT CC (c = a or b). In other words, as long as an orbit of system ()21|) stays near L c , the 
change in h can be estimated by the inequality 

^ > e{A u v c cos6 - 6). (32) 

Recall, that the constants v a and w& are given by (fT2]) . By (fT3]l we may fix the choice of L a and 
Lb such that 

«a > v b . (33) 

Formally, system ([2"T]) is not of the form studied in [7]. However, the study in [7] (see Theorem 
1 there) is reduced to the study of Poincare maps n cc / of the form which includes ([29]) as a partial 
case. Therefore, we may apply the results of [7] to the study of the behaviour of the orbits of 
system (|2"Tj) which belong to the invariant manifolds C^(e) (these orbits intersect the cross-sections 
E Q U so they are defined by the Poincare maps H CC ')- Namely, we have two main conclusions. 
First, we consider the codes £ such that the symbols a and b go always in blocks of some fixed size, 
large enough for the transitions from a neighbourhood of L a to a neighbourhood of Lb or from a 
neighbourhood of Lb to a neighbourhood of L a to happen relatively rarely, so the contribution of 
these transitions to the change in h can be neglected (more precisely, it is absorbed in the small 
6 term in (|34[) ). Then, by (f3"2")l . we obtain that for the orbits that stay on the invariant manifolds 
(e) with such codes £ the evolution of h is estimated by 

dh , , , 
— > e{A u v c{t) cos 9-5) (34) 

where c(t) — a or b indicates where the orbit finds itself at the moment t, near L a or near Lb- 
Second, by repeating the corresponding construction in [7j, for each small e we can find a code £ 
and an initial condition on the corresponding manifold (e) such that for the corresponding orbit 

c(t) = a when A UJ cos9>0, and c(t) = b when A LJ cos9<0. 

For such orbit, we can rewrite (1341) as 



— > e(A UJ cos9 — h lA^cos^l ^ 

do Z Z 

By integrating over one period of 9 = ujt, we find that for this particular orbit the change of the 
energy h over each consecutive 2tt/uj interval of time is estimated by 

Ah>e(\A ul \(v a -v b )-2Tr6)>0 

(recall that 6 can be taken as small as we need and v a > Vb by assumption of the theorem). Thus, 
the energy of the chaotic oscillator steadily grows along the chosen orbit (at a linear rate), which 
completes the proof of the theorem. 



9 



Acknowledgements 

We are grateful to V.Gelfreich, A.Vladimirov, and V.Rom-Kedar for useful discussions. This work 
was supported by the Levcrhulmc Trust grant RPG-279 and the Russian Ministry of Education 
and Science grants 14.B37.21.0862. 

References 

[1] A. Komech, H. Spohn, M. Kunze. Long-Time Asymptotics for a Classical Particle Interacting 
with a Scalar Wave Field. Communications in Partial Differential Equations 22 (1997) no. 1/2 
307-335. 

[2] V. Imaikin, A. Komech, H. Spohn. Scattering Theory for a Particle Coupled to a Scalar Field. 
DCDS-A. 10 (2004) no. 1/2 pp 387-398. 

[3] D. V. Treschev, Oscillator and thermostat, DCDS-A, 28:4 (2010), 16931712 

[4] S. Bolotin, D. Treschev. Unbounded Growth of Energy in Nonautonomous Hamiltonian Sys- 
tems. Nonlinearity 12 (1999) 365-388. 

[5] A. Dclshams, R. de la Llave, T.M. Seara, Ageometric approach to the existence of orbits with 
unbounded energy in generic periodic perturbations by a potential generic geodesic flows on 
T2, Commun. Math. Phys. 209, 353392 (2000) 

[6] A. Delshams, R. de la Llave, T.M. Seara, Orbits of unbounded energy in quasi-periodic per- 
turbations of geodesic flows. Adv. Math. 202, 64188 (2006) 

[7] V. Gclfrcich, D. Turaev. Unbounded Energy Growth in Hamiltonian Systems with a Slowly 
Varying Parameter. Commun. Math. Phys. 283 (2008) pp 769-794. 

[8] V.I. Arnold, V.V. Kozlov, A.I. Neishtadt, Mathematical Aspects of Classical and Celestial 
Mechanics. Encyclopaedia of Mathematical Sciences, vol. 3. Springer Verlag, 1988. 

[9] V.I. Arnold, Mathematical methods of classical mechanics. New York: Springer- Verlag, 1989. 

[10] L.P. Shilnikov, On a Poincarc-Birkhoff problem, Math. USSR Sb. 3, 91102 (1967) 

[11] V.S. Afraimovich, L.P. Shilnikov, On critical sets of Morse-Smale systems, Trans. Moscow 
Math. Soc. 28, 179212 (1973) 

[12] L.P. Shilnikov, A.L. Shilnikov, D. Turaev, L.O. Chua. Methods of Qualitative Theory in 
Nonlinear Dynamics, Part I. World Scientific, 2001. 

[13] L.P. Shilnikov, A.L. Shilnikov, D. Turaev, L.O. Chua. Methods of Qualitative Theory in 
Nonlinear Dynamics, Part II. World Scientific, 2001. 

[14] D. Turaev, V. Rom-Kedar. Elliptic Islands Appearing in Ncar-Ergodic Flows. Nonlinearity 
11 (1998) 575. 

[15] Ya. G. Sinai. Dynamical Systems with Elastic Reflections. Ergodic Properties of Dispersing 
Billiards. Russ. Math. Surv. 25 (1970) pp 137-189. 

[16] L. A. Bunimovich, Ya. G. Sinai, N. I. Chernov, Markov partitions for two-dimensional hyper- 
bolic billiards, Russian Mathematical Surveys, 1990, 45:3, 105152. 



10 



[17] A. Katok, J. Strelcyn. Invariant Manifolds, Entropy and Billiards; Smooth Maps with Singu- 
larities. Springer- Verlag, 1986. 

[18] N. Chernov, R. Markarian. Chaotic Billiards. American Mathematics Society, 2006. 

[19] M.W. Hirsch, S. Smale, R.L. Dcvancy. Differential Equations, Dynamical Systems & An 
Introduction to Chaos (2 nd Edition). Elsevier Academic Press, 2007. 

[20] N. Fenichel, Persistence and smoothness of invariant manifolds for flows, Indiana Univ. Math. 
J. 21, 193226 (1971). 



11 



